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I.  INTRODUCTION 


1  he  parameter  space  method  provides  an  analytical  tool  developed  for  use  in  control  system 
analysis  and  synthesis  Although  not  necessary,  its  application  is  facilitated  by  augmenting  the 
analytical  results  with  graphical  portrayals  in  a  selected  multi-parameter  space.  The  method 
requires  that  the  control  system  be  described  by  a  characteristic  equation  which,  for  sarnpled- 
data  or  digital  systems,  may  be  expressed  in  the  /-domain.  The  technique  is  based  on  the 
analysis  and  synthesis  methods  for  linear  and  nonlinear  control  system  design  which  are 
amply  described  in  Siljak’s  excellent  monograph  on  the  subject.'  In  another  work,  Siljak 
describes  the  application  of  the  technique  to  the  analysis  and  synthesis  of  linear  sampled-data 
control  systems. : 

Once  the  system  characteristic  equation  has  been  obtained,  the  parameter  plane  method 
enables  the  designer  to  evaluate  graphically  the  locations  of  roots  of  the  equation.  Hence,  he 
may  design  the  control  system  in  terms  of  the  chosen  performance  criteria;  e.g.  absolute 
stability,  damping  ratio,  and  setting  time.  He  is  able  to  see  the  effect  on  the  characteristic 
equation  roots  of  changing  two  adjustable  parameters.  Siljak  further  simplified  the  design  by 
introducing  Chebyshev  functions  into  the  equations,  thereby  putting  them  in  a  form 
that  is  amenable  to  their  solution  by  a  digital  computer. 

The  method  has  been  extended  to  portray  the  effect  of  varying  the  sampling  period.1 2 3’4  The 
extended  method  permits  one  to  see  the  effect  of  the  choice  of  values  assigned  to  the  sampling 
period  on  absolute  and  relative  stability.  Also,  the  recursive  formulas  shown  therein  are 
simpler  in  form  than  the  Chebyshev  functions  of  Reference  2.  The  resulting  formulation  is 
deliberately  cast  in  a  form  that  makes  it  particularly  amenable  to  solution  by  a  digital 
computer  or  a  desk  calculator,  again  emphasizing  the  interplay  between  analysis  and 
computing  machines.  When  portrayed  graphically,  the  results  show  the  dynamic  relation 
between  the  selected  parameters  and  the  characteristic  equation  roots,  as  a  function  of  the 
nondimensional  independent  argument,  <unT.  Hence  one  readily  can  deduce  the  dynamic 
effect  upon  the  system  of  various  combinations  of  values  of  the  selected  parameters  defining 
the  parameter  space. 


1  II  I).  Siljak.  Nonlinear  Systems.  Wiley.  New  York.  1969 

2  D.D.  Siljak.  Analysis  and  Synthesis  of  feedback  Control  Systems  in  the  Parameter  Plane.  Part  ll-Sampled  Data  Systems." 
Transactions  of  the  Institute  of  electrical  anti  f  leet ronics  Engineers.  Part  II  Applications  and  Industry,  Vol  83,  November 
1964,  pp,  458-466 

3.  S.M  Selt/er.  "Sampled-Dala  Control  System  Design  in  the  Parameter  Plane."  Proceedings  of  the  Ath  Annual  Allerton 
Conference  Circuits  and  System  Theory.  1970.  pp.  454-461 

4  S  M  Selt/er.  "fnhancing  Simulation  Ffficicncy  with  Analytical  Tools,  ~  Computers  and  Electrical  Engineering,  Vol.  2, 1975. 
pp.  35-44. 
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I  he  history  of  the  continuous-lime  domain  version  of  the  parameter  plane  technique  is  well- 
descrihed  with  suitable  references  in  Reference  I  Briefly  summarizing  that  history,  in  1876 
I  A.  Vishnegradsky  of  the  Leningrad  School  of  Fheoretical  and  Applied  Mechanics  developed 
and  used  the  First  version  of  the  parameter  plane  technique  to  portray  system  stability  and 
transient  characteristics  of  a  third  order  system  on  a  two-parameter  plane.  In  1949  Professor 
Yu  l.  Ncimark  of  the  Russian  School  of  Automatic  Control  generalized  Vishnegradsky 's 
approach  to  permit  the  decomposition  of  a  two-parameter  domain  (D)  describing  an  nth  order 
system  into  stable  and  unstable  regions  The  technique  was  called  ^-decomposition.  During 
the  period  1959-1966  Professor  D  Mitrovic.  founder  of  a  Belgrade  group  of  automatic 
control,  extended  the  method  to  enable  the  analyst  to  relate  the  system’s  variable  parameters 
to  the  system  response,  using  the  last  two  coefficients  of  an  nth  order  characteristic  equation. 
Beginning  in  1964.  Professor  D  D.  Siljak.  then  a  student  of  Mitrovic's  at  the  University  of 
Belgrade,  generalized  the  method  and  called  it  the  Parameter  Plane  method.  His  method 
permitted  the  analyst  to  select  an  arbitrary  pair  of  characteristic  equation  coefficients  (or 
parameters  appearing  within  the  coefficients)  and  portray  both  graphically  and  analytically 
the  dependence  of  the  system  response  upon  the  selected  parameters.  The  method  was 
extended  subsequently  by  Siljak  and  others  to  encompass  a  host  of  related  problems.  In  1967 
Professor  J.  George  modified  the  D-dccomposition  method  to  enable  the  portrayal  of  the 
absolute  stability  region  in  a  multi-parameter  space.  George  also  showed  how  to  portray 
contours  of  relative  stability,  as  did  Siljak.  All  of  the  foregoing  work  is  carefully  and 
completely  referenced  within  Reference  I.  In  1966  and  subsequent  years.  Seltzer  has  applied 
the  parameter  space  method  to  the  design  of  missile,  aircraft,  and  satellite  controllers, 
including  systems  containing  one  or  two  nonlincarities;  the  analysis  of  the  dynamic  effects  of 
the  uonliner  “Solid  Friction"  (Dahl)  model  for  systems  with  ball  bearings,  such  as  control 
moment  gyroscopes  and  reaction  wheels;  and  the  specification  by  the  system  designer  of  the 
dynamic  structural  flexibility  constraints  to  the  structural  designer.  Most  of  the  work  has 
appeared  in  the  technical  journals  of  the  Institute  of  Electrical  and  Electronic  Engineers 
(IEEE),  the  American  Institute  of  Astronautics  and  Aeronautics  (AIAA),  the  International 
Journal  of  Control,  and  the  journal.  Computers  &  Electrical  Engineering.  A  portion  of  the 
history  that  has  not  been  reported  upon  previously  (with  one  exception  to  be  noted)  is  the 
control  system  work  conducted  by  the  German  rocket  scientists  in  the  early  I940’s  at 
Peenemunde.  There,  Dr.  W.  Haeussermann  and  others  applied  the  D-decomposition 
technique  to  the  design  of  the  V-2  Rocket,  following  Dr.  Haeussermann’s  earlier  (prc-World 
War  1 1 )  application  to  the  control  of  an  underwater  torpedo.  This  work  was  not  published  in 
the  open  literature  because  of  national  security  constraints.  When  the  group  came  to  the 
United  States  to  work  with  the  Army  Ballistic  Missile  Agency  (first  at  Ft.  Bliss,  Texas,  then  at 
Redstone  Arsenal.  Huntsville,  Alabama).  Dr.  Haeussermann  and  his  associates  continued  to 


apply  the  method  to  US  Army  missiles  (and  later  to  NASA  space  vehicles).  Again,  national 
(this  time,  another  nation!)  security  precluded  publication  in  the  open  literature  until  1957.5 

In  1964  Professor  Siljak  published  the  first  application  of  the  parameter  plane  technique  to 
sampled-data  systems.  As  mentioned  abo\e.  this  was  extended  in  References  3  and  4.  In  1971 
Seltzer  presented  an  algorithm  for  systematically  solving  the  Popov  Criterion  applied  to 
sampled-data  systems  '1  Applications  of  these  sampled-data  parameter  space  techniques  are 
found  in  References  3,  4  and  6. 

2.  ANALYTICAL  PRELIMINARIES 

The  technique  requires  that  the  control  system  be  described  by  a  characteristic  equation  in 
the  z-domain.  Two  adjustable  parameters  ( ko.  k] )  are  selected,  and  the  characteristic  equation 
(CE)  is  recast  in  terms  of  them:  i.e. 


CE  -  Z  Y  ;cj  «  0, 
j-0  J 

(1) 

"  Yj  (Jj’  Ro’  V 

(2) 

Ts  16 

(3) 

z  *  e  =  re  , 

r  -  r(S,con,T)  =  e  ^WnT’ 

(4) 

B  =  B(£,u)  ,T)  =  cos0=  cos  (uj  T  v/l =  cosuiT 
n  n 

(5) 

where  £,  a*,  and  T  represent  the  damping  ratio,  natural  frequency,  and  sampling  period, 
respectively.  The  symbol  d,  represents  all  system  parameters  other  than  k„,  ki.  To  transform 
the  characteristic  equation  from  an  nth  order  polynomial  into  an  algebraic  equation,  z’  may  be 
defined  in  terms  of  its  real  and  imaginary  parts,  R,  and  I,: 

zJ  =  Rj  +  1  li  ‘  (6) 

5.  W.  Haeussermann.  “Stability  Areas  ol  Missile  Control  Systems."  Jet  Propulsion.  July  1957.  pp.  787-795. 

6.  S.M.  SelOer.  “An  Algorithm  lor  Solving  the  Popos  Criterion  Applied  to  Sampled-I)ata  Systems."  Proceedings  of  the  Third 
Southeastern  Symposium  on  System  Theory.  5-6  April  1971.  pp.  G5-0  through  G5-7. 
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It  readily  follows  In  in  I  quaiion  (J)  that  the  values  ol  K  and  I  are 

Rj  *  rJ  cos  JO  (7R) 

and 


I 


J 


sin  JO  . 


(71) 


While  equations  ( ’’R  t  and  ( 7 1 1  arc  satisfactory  lor  determining  the  values  of  R,  and  l„  a  set  of 
recursive  formulas  can  he  derived  that  arc  particularly  amenable  to  implementation  on  a 
desktop  calculator  or  digital  computer.  They  are  shown  below  as  1  quation(K  land  are  derived 
in  the  appendix 


J  +  l 


-  2rBX 


J 


r2X 


J-l 


(8) 


where  X  may  be  either  R  or  l,.  Only  two  values  of  \,  (i  e  .  two  values  each  of  R,  and  I,)  are 
needed  to  obtain  iteratively,  one  at  a  time,  all  other  values  ol  R  .  I  I  hesc  are  obtained  from 
the  definition  ol  /  in  Fquation  (3)  When  |— 0.  the  value  of  t  is  unity  and,  from  Fquation  (b). 


Z°- 


R,.  +  l  I, 


(9) 


or 


1  -  R 
0  -  I 


0 

0  * 


When  j— I,  the  value  of  /  is.  in  Fulcr  form. 
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z  «  r  cos  0  +ir  sin  0 


R1  +  1  L1 


Hence. 
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( 10R) 

(101) 


(11) 
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W. 


rB 


(12R) 


a  mi 
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1 


(121) 


ll  will  turn  out  to  be  useful  in  the  sequel  to  obscr\c  from  Equation(K)  that  the  radical,  V  I^B\ 
appears  as  a  factor  in  each  \alue  of  I  Solely  lor  simplicity,  it  may  he  factored  out  hy  defining  a 
new  term.  I  as 


I 


i 


V 


(13) 


It  may  he  observed  that  K  and  I  are  functions  of  r  and  B  or  of  £.cu„,  and  1  1 1  Equation  (6)  is 
substituted  into  Equation  ( 1 1  and  the  real  and  imaginary  parts  of  the  resulting  equation  are 
separated,  two  simultaneous  algebraic  equations  are  obtained 

These  two  equations  contain  the  adjustable  parameters,  or  variables.  It,,.  It,.  Hence  it  maybe 
observed  that  the  two  equations  may  be  solved  explicitly  for  k<>  and  k,  as  functions  of  the 
other  system  parameters  (d,t  and.  in  particular,  as  functions  of  the  independent  argument. 
<o.,  I .  It  is  the  latter  observation  that  forms  the  basis  ol  the  parameter  space  method  for 
determining  stability  regions  and  system  dynamic  response  in  terms  of  selected  system 
parameters. 


.V  STABILITY  DETERMINATION 


I  he  method  involves  the  definition  of  stability  boundaries  on  a  multi-parameter  (to  include 
k„.  k  i )  space.  I  hesc  boundaries  arc  found  from  the  pair  of  simultaneous  algebraic  equations 
that  result  from  the  following  operations  on  the  system  characteristic  equation  written  in 
terms  of  the  complex  variable,  /.  Equation  (b)  is  substituted  into  characteristic  Equation  ( 1 ). 
and  the  resulting  real  and  imaginary  parts  are  separated  and  equated  to  zero.  The  two  resulting 
algebraic  equations  may  be  solved  for  any  two  parameters  (such  as  ko  and  ki)  within  them. 


There  may  be  as  many  as  four  stability  boundaries  for  any  system  described  by 
characteristic  Equation  (I),  although  it  is  not  necessary  for  all  four  to  exist. 

One  stability  boundary  separates  the  stable  complex  conjugate  pairs  of  roots  from  the 
unstable  ones.  It  consists  ol  a  map  of  the  unit  circle  from  the  complex  /-plane  onto  a  selected 
(such  as  ko.  kit  parameter  space.  Initially  one  may  consider  the  oft-occurring  case  where  the 
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coefficients  >,  of  the  powers  of  /  in  the  characteristic  equation  are  linear  combinations  of  k«, 
Wi.  i.e.. 


YJ  * 


k  +  b.  k.  +  c 
o  j  1 


(14) 


where  a„  b„  c,  represent  all  system  parameters  other  than  ko,  ki.  In  this  case  the  two 
simultaneous  equations  resulting  from  the  real  and  imaginary,  respectively,  parts  of  the 
characteristic  equation  assume  the  form. 


Re  |C.E.|  -  Aj  a  +  b-C^  -  0, 

(15R) 

Im  |C.E.|  -  a  +  b-Cj  " 

(151) 

ki  readily  by  applying  Cramer’s  Rule,  yielding 

(16) 


and  may  be  solved  for  ko  and 


B1C2  "  B2C1 


A2C1- 


A1C2 


(17) 


J  -  A1B2  '  A2 

B1 

(18) 

n 

n 

n 

A1  “  Z  aiR1 
j»o  3  J, 

B1 

■ z  bjRi 

j-0  J  \ 

C  -  I  c  R 
j-0  J  \ 

n 

n 

n 

A2  "  1  Vj 
1-0  J 

B2 

'  E  bixi 
3-0  J 

C2  ‘  E  C1X1 

L  j-o  J  \ 

(19) 

where  J  represents  the  Jacobian  associated  with  the  pair  of  algebraic  equations.  To  find  the 
stability  boundary  in  question,  one  merely  sets  £  equal  to  zero  in  the  definitions  of  R,  and  I, 
used  in  obtaining  k<,  and  ki  in  Equations  ( 16)  and  ( 17).  The  result  is  a  boundary  in  the  k0-ki 
parameter  plane  that  is  defined  in  terms  of  system  parameters  dj  [in  the  general  case,  see 
Equation  (2)]  or  a,.  b„  c,  in  the  special  case  of  Equation  (14)  and  the  independent  argument, 
<u„T.  The  independent  argument  is  varied  in  value  between  zero  and  n  (assuming  that  the 
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system  being  designed  possesses  low  -pass  tiller  characteristics  so  that  only  the"primary  strip” 
need  be  considered  I.  thereby  defining  the  boundary  in  terms  ol  k,  and  k  It  the  coefficients  a  , 
b.  and  c,  contain  exponential  terms  with  I  appearing  in  the  exponents,  then  each  exponent 
may  be  replaced  by  its  power  series  and  truncated  according  to  the  accuracy  that  is  desired  If 
the  two  pa  i  a  meters  k  and  k  do  not  appear  linearly  as  expressed  in  Equation  (141.  they  may  be 
solxed  tor  h\  substituting  Equation  (6)  into  the  characteristic  Equation  ( 1 1  and  separating  the 
result  into  the  real  and  imaginary  parts  1  his  still  results  in  two  simultaneous  algebraic 
equations  which  may  be  solved  for  the  two  variables,  k,  and  ki,  although  not  as  handily  as 
expressed  in  Equations  (lb)  and  ( 17) 

I  he  next  two  stability  boundaries  are  those  separating  the  stable  real  roots  from  the 
unstable  ones  These  comprise  a  mapping  ol  the  /  I  and  z=-l  points  from  the  /-plane  onto  the 
selected  parameter  space.  These  are  found  readily  by  substituting  /=  1  and  /=-l ,  respectively, 
into  the  characteristic  Equation  1 1 )  I  ach  ol  the  two  resulting  equations  results  in  a  definition 
ol  the  two  real  root  boundaries  in  the  selected  parameter  space. 

The  fourth  boundary  is  a  mapping  of  the  conditions  that  cause  the  two  simultaneous 
algebraic  equations  (used  to  define  the  complex  conjugate  root  boundary)  to  become 
dependent  In  the  case  where  k<>  and  k:  appear  linearly  as  in  Equation!  14).  this  case  is  found  by 
determining  the  conditions  that  cause  the  Jacobian  of  Equation  (IS)  to  become  identically 
equal  to  zero 

For  a  linear  samplcd-data  control  system  to  be  stable,  it  is  necessary  that  all  roots  of  the 
characteristic  equation  lie  within  (he  unit  circle  on  the  /-plane  It  the  system  is  low-pass  in 
nature,  the  stable  region  is  bounded  by  the  semi-circle  defined  by  the  upper  half  of  the  unit 
circle  (as  discussed  above)  and  the  singularities  associated  with  /  =  ±1  (as  discussed  above) 
and  .1=0  (as  discussed  above)  T  he  mapping  of  these  boundaries  onto  the  parameter  space  will 
bound  the  stable  region,  if  one  exists.  1  hat  region  is  determined  by  applying  a  shading 
criterion  or  using  a  test  point  If  the  Jacobian  is  greater  than  zero,  then  the  stable  region  (if  it 
exists)  lies  to  the  left  of  the  complex  conjugate  root  boundary  as  i«„T  increases;  the  left  side  of 
the  line  is  double  cross-hatched  to  indicate  a  boundary  associated  with  double,  or  complex 
conjugate  roots.  (If  the  Jacobian  is  less  than  zero,  the  stable  region  lies  to  the  right).  Single 
cross-hatching  is  used  on  the  two  contours  associated  with  the  real  roots.  The  side  of  the 
contour  on  which  to  place  the  cross-hatching  is  determined  by  the  requirement  that  cross- 
hatching  be  continuous,  or  on  the  same  side,  of  the  contours  as  the  intersections 
corresponding  to  /=  + 1  and  /=  -I  arc  approached  along  either  the  complex  root  or  real  root 
stability  boundary.  In  the  unusual  (but  physically  possible)  singular  case  when  .1=0.  the 
shading  of  the  resulting  boundary  is  determined  in  a  similar  manner,  i.e..  the  side  on  w  hich  to 


shade  the  boundary  in  determined  by  the  physical  requirement  that  the  number  m  stable  roots 
must  nescr  become  less  than  zero  as  one  moves  across  a  stab.  I  us  boundary  1  his  J=0  case  may 
arise  lor  a  particulai  frequency ,  o»,  or  lor  a  particular  combination  ol  system  pa.  ameters  The 
latter  situation  is  included  in  the  work  concerning  a  spinning  Sky  lab 

4.  SYSTEM  DYNAMICS 

Once  the  stable  region,  it  it  exists,  has  been  determined  in  the  selected  parameter  space,  the 
dynamic  or  transient  characteristics  of  the  sy  stem  can  be  specified  in  terms  of  the  locations  of 
the  roots  of  the  characteristic  equation  (pole  placement).  For  the  complex  conjugate  roots, 
these  locations  are  defined  in  terms  of  damping  ratio  (£)  and  system  natural  frequency  (w„). 
Contours  of  constant  £  are  determined  as  functions  of  mn  and  T  in  precisely  the  same  manner 
that  the  complex  conjugate  stability  boundary  was  determined  except  that  £ is  not  set  equal  to 
zero. 


The  real  root  locations  corresponding  to  values  of  /  when  6  equals  0°  and  180°  also  may  be 
plotted  on  the  parameter  plane  by  setting  /  equal  to  a  positive  or  negative  real  constant, 
substituting  that  value  into  Equation  ( I ).  and  solving  for  k>  as  a  function  ofko.  Each  resulting 
contour  corresponds  to  a  location  of  a  real  root  in  the  z-domain.  When  z=+a 


l  Y,aJ  -  0.  (20) 

J-0  J 


and  when  z=  -a. 


n/2 

t  y 
J-o 


2) 


2  j 
a 


(n-2)/2  (2j+li) 

E  Y(2.+1.)oi  =  0,  n  even 

j-0  J  J 


(n 


D/2 
l  (Y 
j-0 


a2j_  a(2j+D, 

2j“  T(2j+n 


0,  n  odd 


(2D 


(22) 


where  a  is  a  positive  real  number  (The  two  associated  stability  boundaries  may  be  found  by 
setting  a  equal  to  unity).  Y'alues  of  o  corresponding  to  desired  real  root  locations  may  be 
substituted  into  Equations  (20)  through  (22)  to  obtain  parameter  space  contours 
corresponding  to  these  locations. 


7  S  M  Seltzer.  “Passive  Stability  of  a  Spinning  Sky\*b.m  Journal  of  Spacecraft  ami  Roc  krts,  Vol.9,  No  9.  September  W*T2.pp 
651*655 
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Now  the  dynamic  ef  fect  of  any  chosen  design  point  in  the  parameter  space  may  be  specified 
in  terms  of  £,  <u„,  a,  and  T.  1  he  analytical  technique  developed  permits  the  designer  to  observe 
the  effect  of  simultaneously  changing  three  control  parameters  and  the  sampling  period.  Most 
existing  conventional  techniques  permit  the  observation  of  the  effect  of  changing  only  one 
control  parameter  and  don’t  show  the  effect  of  various  sampling  periods.  If  the  designer  is 
clever,  sometimes  the  system  parameters  may  be  manipulated  within  the  equations  defining 
the  stability  and  dynamic  contours  so  that  more  than  two  parameters  (such  as  ko,  ki)  may  be 
used  to  define  a  parameter  space.  An  example  of  a  three-parameter  space  is  provided  in  the 
sequel. 

Sometimes  it  is  specified  that  the  setting  time  of  the  system  be  less  than  a  prescribed  value. 
This  corresponds  to  requiring  that  the  real  part  of  the  roots  of  the  characteristic  equation  be 
less  than  a  prescribed  negative  real  constant.  A  boundary  corresponding  to  this  requirement 
can  be  drawn  on  the  parameter  plane  by  mapping  a  circle  of  constant  radius  (for  a  chosen 
values  of  <o„,  £,  and  T)  from  the  z-plane  onto  the  parameter  plane.7  Relations  exist  for 
estimating  the  maximum  overshoot  and  peak  time  of  transient  response  when  it  is  valid  to 
assume  a  second  order  system."  However,  a  simple  estimate  sometimes  can  be  made  merely  by 
looking  at  the  difference  equation  representing  the  system  response,  estimating  when  the 
overshoot  will  occur,  and  plotting  a  corresponding  line  on  the  parameter  plane.  Steady  state 
responses  may  be  found  from  the  open  loop  transfer  function  and  the  assumed  forcing 
functions  in  the  conventional  manner. 

5.  EXAMPLE 

Consider  the  planar  model  that  portrays  controlled  rigid  body  rotational  dynamics.  It  is 
shown  in  block  diagram  form  in  Figure  1.  The  plant  to  be  controlled  is  represented  by  the 
transfer  functions, 

G4(s)  =  <&(s)/Tc(s)  =  l/dvs  (23) 

and 

G5(s)  =  $(s)/<S>(s)  =  1/s  (24) 


8.  B.C.  Kuo.  Analysis  and  Synthesis  of  Sampled- Data  Control  Systems.  Prentice-Hall.  New  Jersey,  1963. 
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Figure  1.  Block  diagram  of  controlled  rigid  body  motion. 

where  s  represent;*  the  l  aplaee  operator,  </>  represents  the  iigtd  both  at t mule  [in  the  time 
tlomnin.  </Ht),  in  the  I  aplaee  donum.  dH st|,  1,  represents  the  eommundetl  totquc,  and  the 
overdot  represents  the  derivative  with  respect  to  time.  !  he  two  sensors  that  measure  altitude 
td>'  and  attitude  rates  (d>f  arc  assumed  to  be  perfect  with  unity  transfer  functions  1  he  on- 
boaid  digital  controllei  develops  the  commanded  control  torque  I rom  the  input  states  d*  and  d> 
and  the  commanded  attitude  signal  </>..  tih„(sl  represents  a  /ero-ordei  hold  in  (lie  computer 


Ch  («)  -  (l-o  'S)/s, 
o 


where  1  is  the  sample  period  of  the  on  board  digital  computet  Hie  l'll>  control  algorithm  is 
represented  by  (iits)  and  Ch(sl 

C.jCh)  -  Kr  +  K,/m  (2b) 


where  Ki>aiul  h.i  are  the  position  and  integial  teed  back  gains  and 


whore  K\.  in  i ho  ilciivativc  Icoilback  )Mmi  imoio  coinmonlv  i.illcil  the  into  |imnl  I  hr  llnul 
oulct  chaiacteiiNtn  equation  in  I  quation  i  I i  with  n  >  where  the  corllicienlN  ol  /'  air  in  I  hr 
loim  ol  I  ijnalion  i  I4i.  io  . 


>0  -  (a  -  h  t  -  n. 
>.  -  (-2.1  f  2c  4  \), 

-  (a  4  b  O  . 

£ 

>3  -  l. 


Modified  gams  a.  b,  c  me  defined  as 


krt/jv 


b  a  K  TV2J  , 
P  v 

c  »  K.T'V^J  , 

v 


( 2  H) 
(2^) 
(  H>) 
(  U> 

('2) 
(II) 
(  O) 


I  lie  m abilns  botimlai  ion  are  piesentcd  in  lot  iiin  ol  a.b.o.  ami  <nn  I  In  obsei ving  llial  the  t  ootN 
(m  iIiin  ease  tluee)  ol  a  characteristic  o»|iiaiion  irpirsenting  a  Niablo  linoai  Namplovl  data 
oontiol  N\Ntom  nuist  all  ho  within  the  unit  circle  in  the  /  plane  It  in  assumed  that  the  nvnIciii 
poNNCNNCN  Ion  paNN  liltci  chaiacleristics  so  that  onlv  thcpmuniv  stnptcoiicspondmgtoO  •  I* 
•  rr)  need  be  considered  I  lie  stable  region  niav  ho  defined  In  mapping  Kn  tluee  bounilaiiCN 
horn  the  >  plane  onto  the  a.  b.  c  paimnctci  space  b\  liiNt  consntoiing  paiaineteis  a.  bto  be  k... 
ki  in  I  ^nation  1 14 1  I  he  hotnulatv  at  /  H  in  louinl  b\  substituting  that  value  tot  /  into  the 
('  I  .  Molding  the  Ntnbililv  bomularv 

c  *  0  (  VO 

I  he  bomularv  at  i  I  in  IoiiihI  b\  substituting  that  value  lot  ;  into  the  t'  I  .  Molding 

a  -  2 .  0*0 


I  he  complex  couitigafe  toot  stabilitv  boumlaiv  mav  be  lomnl  In  setting  /  e»|iial  to  cos  ,ii„  I 
*  i  sin  in,,  I  in  the  C  I  llowovoi.  at  this  point  it  is  convenient  to  use  the  iociumvo  (onnulas  ot 
I  ipiation  (St  ami  tiaiistouu  the  t '  I  into  two  alpebi.no  cvpiations  bv  scpataiiug  the  teal  ami 


imaginary  parts.  If  they  are  solved  lor  a  and  hand  the  associated  Jacobian  J,  one  obtains  as  the 
complex  conjugate  root  stability  boundary. 


a  =  CO  +  B)  /  ( 1  -  B)]  c  +  (1  -  B),  (37) 

b  -  c  +  (1  -  B),  (38) 


and 


J  «  -4  (1  -  B)  <  0  V  0  <  0  <  it. 


(39) 


Examination  of  Equation  (39)  reveals  that  the  stability  boundary  associated  with  the  singular 
case  J=0  only  occurs  when  B-l,  which  is  the  already  considered  /=  I  case.  Now  the  stability 
boundaries  of  Equations  (35).  (36),  (37)  and  (38)  can  be  plotted  on  a  three-dimensional  plot 
with  axes  a,  b,  c.  A  sketch  of  these  boundaries  is  shown  in  Figure  2.  The  stability  region  is 
found  by  applying  the  “cross-hatching"  rules  or  by  using  one  or  more  test  points(known  to  be 
stable,  to  lie  on  a  stability  boundary,  or  to  be  unstable)  1 1. IX)  the  stable  region  (if  it  exists)  lies 
to  the  left  of  the  complex  conjugate  root  stability  boundary  of  Equations  (37)  and  (38),  as  6 
(or  tu„T)  increases;  the  left  side  of  the  line  is  double  cross-hatched  to  indicate  a  boundary 


{a  -  b  -  2c  [1  -  (b-c ) -1 )  } 


Figure  2.  Stability  boundaries  in  3-dimensions. 


associated  with  double.  01  complex  conjugate,  roots  II  l<().  as  in  this  case,  the  stable  region 
lies  to  the  right  and  the  tight  side  ol  the  line  is  double  cross-hatched  Single  cross-hatching  is 
used  on  the  contours  ol  I  quations  (35)  and  (3<u  associated  with  single  root  boundaries.  The 
side  ol  the  contour  on  w  Inch  to  place  the  ct  oxs-hatching  is  deter  mined  b\  the  requirement  that 
cross-hatching  be  continuous,  or  on  the  same  side,  ol  the  contours  as  the  intersections 
corresponding  to  /  t  I  (()  0)  and  /  -ltd  n  l  ate  approached  along  either  the  complex  root 
or  real  root  stability  boundary  losing  this  cr itenon.  the  stable  region  is  lound  to  be  bounded 
by  thee  Oplane.thea  2  plane. and  thecuivespecilied  by  Fquations(37)and(3K).  The  latter 
may  be  solved  lor  a  as  a  function  ol  b.  resulting  in  the  simple  expression 

a-b  -2c  Cl-(b-c)  1  3  .  (AO) 

I  he  limiting  values  lor  a  and  b  may  be  lound  by  letting 0  approach!)  and  rr.  i.c..  /=±|.  For  the 
case  of  6  —  0.  represent  B  by  the  series  expansion 

B  -  cos  0  -  l-(O2/2!)+(04/4!)-(O6/6!)+...  (41) 

Then  using  Equations  (37).  ( 38).  and  (41).  one  obtains 

lim  b  ■  c 

0  >  0  (42) 

11m  a  .  -  lim  c  [(4/02)-(2/3)  +  (O2/36)- _ ]  -*•  <* 

0  -♦  0  0  -*■  0  (43) 

I' or 'c  0.  fqu.it  ion  (42)  still  holds  (recognize  ki  0  is  really  the  case  one  is  interested  in).  1  hen, 
using  Equations  (34)  and  (37).  one  obtains 

lim  a  -  11m  (k  03/4J  ui3)  [2- (02/2  !  )  +  (0A/4  ! ) - ] 

0  -*  0  0  -+  0  v 

(0  /2!)-(0A/4 !)+••. 

-  11m  (kI0/2Jvui3)[2-(2O2/6)+(O4/72) - ]-0 


(A4) 


For  the  limiting  case  ot  6  *  n 


llrn  a  -  11m  {  C(  1+B)c/ (  1-B)3+(J-B)  }  -  2 

0  -*  TI  0  -»  V 


(45) 


and 


lim  b  -  lim  (c+l-B)  ■  c+2 

0  -»  n  0  n  (46) 

The  results  of  the  foregoing  are  sketched  graphically  in  two  dimensions  on  Figure  .1,  for  the 
cases  where  c  =  0  and  when  c  *  0. 

I  he  design  technique  employs  the  pole  placement  approach.  It  is  based  on  the  premise  that 
the  dynamic  behavior  of  the  system  is  related  closely  to  the  location  of  the  roots  of  its 
associated  characteristic  equation.  I  he  method  shows  both  analytically  and  graphically  the 
direct  correlation  between  these  roots  and  the  control  gains  and  sample  period  of  the 
controller  The  design  technique  then  involves  the  specification  of  the  characteristic  equation 
root  locations  and  the  subsequent  determination  of  the  control  system  gams  and  sample 
period  needed  to  attain  these  locations.  The  control  system  designer  then  must  determine  the 
system  response  resulting  from  using  these  numerical  values  and  assess  its  adequacy.  If  it  is  not 
adequate,  he  usually  relies  on  his  experience  to  relocate  the  roots  to  improve  the  response  in 
the  manner  desired  for  his  particular  system  (i.e.,  faster  settling  time,  lower  peak  overshoot, 
etc.).  It  is  assumed  that  one  wishes  the  pair  of  complex  conjugate  poles  of  the  C.F.  todominate 
the  dynamic  response  of  the  system  This  response  w  ill  then  be  modified  by  locating  the  third 
(real)  root. 

First  consider  the  pair  of  complex  conjugate  roots.  One  may  use  Equations  (lb)  -  (19), 
substituting  a.  b  for  k(>,  ki,  respectively  |also  using  Equations  (28)  through  (3D).  Upon 
examining  the  nature  of  the  equations  fora,  b,  one  sees  that  they  are  functions  of  r,  B.  and  c,  or 
(finally)  functions  of  (,  0,  c.  For  each  value  of  £,  b  versus  a  may  be  determined  analytically  and 
plotted  graphically  as  a  function  of  6  for  a  given  value  ofc.  I  his  may  be  repeated  fora  number 
of  values  of  c.  Now  the  value  of  0  may  be  selected  (for  selected  values  ofc)  that  best  meets 
whatever  design  criteria  one  chooses.  In  effect,  selecting  0  and  c  will  also  choose  the  real  root 
location.  This  may  be  shown  by  designating  the  real  root  of  the  C.E.  as  ak,  where 

zR  =  6  -  e  -V  (47) 


lb 


NOTE:  NUMBERS  IN 
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Figure  3.  Stability  boundaries  in  2-dimensions 


and  ok  represents  the  damping  tactor  One  may  suhstuutc  Equation  (47 1  into  C.H.  ( I )  and 
solve  lor  b,  obtaining 


b  -  [(l-6)/l+6)]a+[2/(l-6)  c-  (l-6)2/(l+6)]  (48) 

One  sees  that  Equation  (48)  yields  a  straight  line  contour  of  b  vs  a  for  a  given  value  of  c  and  a 
given  value  of  S.  1  Inis,  for  a  given  value  of  c,  where  the  A-contout  ol  Equation  (48)  crosscsthc 
{-contour  of  Equations  (37)  and  (38).  one  obtains  a  value  of  8  for  the  given  value  of  e  and  those 
specified  values  of  {  and  5.  From  Equation  ( I ).  c  may  be  determined  as  a  function  of  the 
characteristic  equation  root  locations: 

c  -  (a2  +  82  -  2a  +  1) (1  -  6)/4,  (49) 

where  a  and  0  are  the  real  and  imaginary  parts,  respectively,  of  the  pair  of  complex  conjugate 
roots  of  characteristic  Equation  (1).  For  a  specified  location  (o,/3)  of  the  pair  of  complex 
conjugate  roots,  a  value  for  c  is  established  by  setting  a  value  for  S  (location  of  the  real  root). 

An  additional  constraint  is  imposed  by  Shannon’s  Sampling  Theorem: 

to  <  to  /2,  (50) 

s 

where  w  is  defined  in  Equation  (5)  and  to,/  2n  is  the  sampling  frequency, 

u)  -  2n/T.  (51) 

s 

From  Equations  (50)  and  (51),  this  constraint  may  be  restated  as 

0  <  it  .  (52) 


The  stable  region  portrayed  in  Figures  2  and  3  may  be  shown  for  several  selected  values  of  c 
on  Figure  4  using  Equations  (lb).  H7)  and  (48)  As  c  increases  in  value,  the  stable  region 
shrinks  in  sire  and  the  {  contours  lie  other  to  the  right.  Whcnc --O.b.  thc{=l  l>/2  contour  lies 
entirely  to  the  right  (and  outside)  of  the  stable  region.  Finally  when  cS?2,  the  stable  region 
disappears  entirely  (for  a  geometric  explanation,  see  Figure  2). 


I  he  application  ol  the  foregoing  technique  will  now  he  summarized. 

Based  on  desned  system  response  characteristics,  tentatively  select  locations  for  the  three 
roots  ol  the  characteristic  equation.  Although  the  zeros  of  the  system  characteristic  equation 
also  allect  the  dynamics,  placement  ol  the  poles  will  dominantly  affect  the  dynamics. 

•  A  root  location  criteria  might  he  to  select  numerical  values  for  4  and  <u„. 

•  Determine  computer  design  constraints.  One  might  be  to  select  as  large  a  value  of  T  as 
possible  to  minimize  on-board  memory  capacity  requirements. 

•  I  hen  qualitatively  determine  how  much  integral  gain  is  desired  for  the  types  of 
.'.sturbance  inputs  expected.  I  his  gives  an  indication  of  the  value  of  c  to  use. 

Using  the  value  of  c  tentatively  established  in  the  above  paragraph,  plot  the  stability 
boundaries  in  the  b  vs  a  parameter  plane,  using  Equations  (36)  through  (38).  Then  plot  the 
desired  4-contour  (as  a  function  ol  argument  iu„  I )  on  the  b  vs.  a  parameter  plane,  using 
Equations  1 16)  and  ( 1 7). 

Using  various  values  ol  6,  plot  5-contours  on  the  same  plot,  using  Equation  (48). 

Find  an  intersection  of  a  6-contour  w  hich  gives  a  desired  value  for  u>„T  (and  hence  T,  since 
wn  has  already  been  specified)  and  6.  In  general,  if  6  —  0.  its  effects  on  the  system  dynamics  is 
small  compared  to  the  effect  of  the  pair  of  complex  conjugate  poles  (placed  by  o>„T,  4). 

Check  the  response  of  the  system  using  the  values  of  a.  b.  c,  and  <U„T  (and  hence  KP,  Kt.  Km, 
and  T)  associated  with  the  selected  intersection.  If  it  is  unsatisfactory,  reiterate  the  above 
procedure,  selecting  another  value  of  c. 

If  the  design  procedure  were  merely  to  specify  the  three  root  locations  (such  as  in  terms  of  a, 
/J.  5,  <o„T.  6),  iteration  would  not  be  required.  Unfortunately,  in  design  practice  one  must 
exercise  engineering  judgement  in  selecting  the  three  root  locations,  giving  rise  to  iterative 
design  procedures  such  as  the  one  outlined  above.  However,  another  possibility  does  exist.  If 
the  real  root  is  placed  near  the  origin  of  the  z-plane  (6  —  0).  and  if  the  pair  of  complex 
conjugate  roots  arc  located  near  the  unit  circle  in  the  /-plane  (i.e..  specifying 4,  and  tu„T),  the 
latter  two  roots  will  dominate  the  system  response.  Then  the  already  developed  tools  for 
specifying  the  system  response  (such  as  determining  explicitly  the  maximum  overshoot  and 
the  settling  time)  of  second-order  systems  may  be  applied  handily.  The  foregoing  procedure 
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may  be  amplified  through  the  application  of  realistic  numerical  values  for  the  system 
parameters 

Assume  a  digital  on-board  controller,  and  further  assume  that  it  is  desired  to  have  a 
controller  natural  frequency  (tu„)  of  six  rad  s.  and  that  the  moment  of  inertia  (J.)  about  the 
single  axis  is  4.3  x  I0J  kg.m' .  Assume  that  integral  control  is  desired  to  drive  to  zero  the  effect 
of  constant  input  disturbances.  This  means  a  non-zero  value  of  c  is  desired  If  a  value  of  £  of 
I  / n/2  is  selected,  and  if  the  value  of  6  is  desired  to  be  kept  as  close  to  zero  as  possible,  it  is 
implied  from  Figure  4c  that,  for  6  >  0,  no  intersections  of  the  £ -contour  occur  when  c  >0.2. 
Hence,  a  value  of  c  =  0.1  is  chosen.  The  smallest  value  of  6  whose  corresponding  6-line 
intersects  the  £-contour,  with  a  reasonable  factor  of  safety,  is  6  =  0.4  The  value  ofco„T=  1.2  is 
on  the  ^-contour  at  the  intersection,  yielding  a  value  of  0.20s  for  T.  Corresponding  values  of  a 
and  bare  1.41  and  0  68.  respectively.  I  his  gives  a  set  of  gain  values  of  1.14  x  10’  Nm/rad,  1.67 
x  10  N-m  rad-s.  and  2.35  x  10  Vm-s  rad  for  kP.  Ki.  and  Ki>.  respectively.  If  response  studies 
show  that  a  larger  value  ol  integral  gain  is  deemed  necessary,  a  larger  value  of  c  (and 
consequently  a  smaller  value  ol  £i  would  be  selected. 

6.  CONCLUSIONS 

An  analytical  method  for  portraying  stability  regions  in  a  selected  oarameter  space  has  been 
shown  for  a  digital  system.  The  method  requires  that  the  system  characteristic  equation  be 
available  and  expressed  in  the  complex  z-domain.  It  also  is  possible  to  apply  pole  placement  to 
obtain  desired  dynamic  characteristics  using  this  modified  parameter  space  technique.  The 
advantage  of  the  technique  over  existing  classical  sampled-data  methods  is  that  the  stability 
and  dynamic  response  characteristics  are  expressed  in  terms  of  several  (rather  than  merely 
one)  selected  parameters.  Also,  the  sampling  period.  T.  need  not  be  expressed  numerically 
before  the  design  technique  begins,  giving  the  system  designer  one  more  degree  of  freedom. 
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APPENDIX 

DERIVATION  OF  RECURSIVE  RELATIONS 


f-t#. 


v 


25 


It  is  possible  to  obtain  algebraic  recursive  relations  for  the  real  and  imaginary  parts  of  the 
complex  variable.  /.  when  it  is  raised  to  the  j'"  power  ( j  is  a  possible  integer).  Erom  Equation 
(3).  /  is  defined  as  the  vector  (or  complex  variable). 


z 


Ts 

e 


r  ei(J=  r  cos  0  +i  r  sin  0. 


(Al) 


Using  Equation  (5)  in  Equation  (Al),  i.c.,  letting  B  =  cos  0,  leads  to 


z 


rB  +  i 


(A2) 


The  real  and  imaginary  parts  of  r'  may  be  defined  as  R,  and  I,,  respectively,  i.e.. 


J 


(A3) 


If  i- !,  Equation  (A3)  becomes  identical  with  Equations  (Al)  and  (A2),  leading  to 


Rj  =  rB 


(A4R) 


and 


If  j=0,  then  Equation  (A3)  becomes 


0 

z 


+  i  I, 


(A4I) 


(A5) 


or 


Rq  =  1  (A6R) 

and 

I0  =  0.  (A6I) 
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PRECEDING  FNQE  ULAJMC 

-  — ■  - - - 


one  squares  the  left  and  right  hand  sides  of  Equation  ( A 1 1.  one  may  obtain 


2  2  120  2  , 

z  -  r  e  -  r  (  cos  20  +  1  sin  20) 

2  2  2 
“  r  (  2  cos  0  -  1)  +  12  r  sin  0  cos  6 

2  o 

-  2  r  cos  0  (cos  0+1  sin  0)  -  r 

-  2rBz  -  r2.  (A7) 

One  may  now  multiply  each  term  of  Equation  (A7)  by  z‘"‘  to  yield 


k+1 

z 


2rBzk 


2 

r  z 


k-1 


(A8) 


Substitution  of  Equation  (A3)  into  Equation  (AS),  letting  j  =  k-l,  k,  and  k+l,  leads  to 

\+i  +  1  xk+i  "  2rB<Rk  +  1  V  -  r2(Rk-i+  ilk-0  •  <A9> 


Separately  equating  the  real  and  imaginary  parts,  respectively,  of  Equation  ( A9)  leads  to  the 
two  equations. 


\+l  “  2rBRk  “  r\-l  (A10R) 

and 

Xk+1  “  2rBIk  -  r2  Xk-1  *  (A10I) 


One  may  observe  the  identical  forms  of  Equations  (A  10)  and  write  them  both  as  a  single 
equation. 


Vi  -  2rBxk  -  -Vi  • 


(All) 


where  Xk  may  be  used  to  represent  cither  Rk  or  U.  Of  course  the  dummy  index  k  may  be 
exchanged  with  j,  yielding  Equation  (8)  of  the  text.  If  one  knows  the  two  values  Xi  and  Xk-i  of 
Equation  (All),  one  may  determine  the  third  value  XkM.  Since  two  sets  of  values  of  Xk  are 


known,  i.e.,  tor  k  =  0  and  k  -  I  |sce  Equations  ( A4)  and  (A6)|,  as  many  values  of  Xk  as  are 
needed  for  the  particular  problem  at  hand  can  be  determined  recursively  from  Equation 
(All). 
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